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Abstract 

We investigate the Abelian monopole condensation in finite temperature SU(2) 
and SU(3) pure lattice gauge theories. To this end we introduce a gauge invariant 
disorder parameter built up in terms of the lattice Schrodinger functional. Our 
numerical results show that the disorder parameter is different from zero and Abelian 
monopole condense in the confined phase. On the other hand our numerical data 
suggest that the disorder parameter tends to zero, in the thermodynamic limit, 
when the gauge coupling constant approaches the critical deconfinement value. In 
the case of SU(3) we also compare the different kinds of Abelian monopoles which 
can be defined according to the choice of the Abelian subgroups. 
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I. INTRODUCTION 



The dual superconductivity of the vacuum in gauge theories to explain color confinement 
has been proposed since long time by G. 't Hooft [|IJ and S. Mandelstam [TJ]. These 
authors proposed that the confining vacuum behaves as a coherent state of color magnetic 
monopoles. In other words the confining vacuum is a magnetic (dual) superconductor. 
This fascinating proposal offers a picture of confinement whose physics can be clearly 
extracted. Indeed, the dual Meissner effect causes the formation of chromoelectric flux 
tubes between chromoelectric charges leading to a linear confining potential. 

Following Ref. H let us consider gauge theories without matter fields. In order to 
realize gauge field configurations which describe magnetic monopoles we need a scalar 
Higgs field @]. In the 't Hooft's scheme the role of the scalar field is played by any 
operator which transforms in the adjoint representation of the gauge group. Let X(x) 
be an operator in the adjoint representation, then one fixes the gauge by diagonalizing 
X(x) at each point. This choice does not fix completely the gauge, for it leaves as 
residual invariance group the maximal Abelian (Cartan) subgroup of the gauge group. 
This procedure is known as Abelian projection ||. The world line of the monopoles can 
be identified as the lines where two eigenvalues of the operator X(x) are equal. Thus, 
the dual superconductor idea is realized if these Abelian monopole condense. Due to the 
gauge invariance we expect that the monopole condensation should manifest irrespective 
to the gauge fixing. In other words all the Abelian projections are physically equivalent. 
However, it is conceivable that the dual superconductor scenario could manifest clearly 
with a clever choice of the operator X(x). It is remarkable that, if one adopts the so 
called maximally Abelian projection 0, then it seems that the Abelian projected links 
retain the information relevant to the confinement ||. 

It turns out that the Abelian projection can be implemented on the lattice so 
that one can analyze the dynamics of the Abelian projected gauge fields by means of non 
perturbative numerical simulations. Indeed, the first direct evidence of the dual Abrikosov 
vortex joining two static quark- ant iquark pair has been obtained in lattice simulations of 
gauge theories |6|-f|. In particular in Ref. f| we considered the pure gauge SU(2) lattice 
theory and found evidence of the dual Meissner effect both in the maximally Abelian gauge 
and without gauge fixing. Moreover we showed that the London penetration length is a 
physical gauge invariant quantity. 

An alternative and more direct method to detect the dual superconductivity relies 
upon the very general assumption that the dual superconductivity of the ground state is 



realized if there is condensation of Abelian monopoles. Thus, according to Ref. [T(| it 
suffices to measure a disorder parameter defined as the vacuum expectation value of a non- 
local operator with non zero magnetic charge and non vanishing vacuum expectation value 
in the confined phase. However, in the case of non Abelian gauge theories, the disorder 
parameter is expected to break a non Abelian symmetry, while the dual superconductivity 
is realized by condensation of Abelian monopoles. As we have already argued, the Abelian 
monopole charge can be associated to each operator in the adjoint representation by the so- 
called Abelian projection ||[5[] . Indeed, the authors of Ref. [|Kj introduced on the lattice a 
disorder parameter describing condensation of monopoles within a particular Abelian pro- 
jection. On the other hand, recent results [1 1] show that the Abelian monopoles defined 
through several Abelian projection condense, suggesting that the monopole condensation 
does not depend on the adjoint operator used in the Abelian projection procedure. This is 



1 



in accordance with the theoretical expectation that monopole condensation should occur 
irrespective of the gauge fixing procedure. However, a gauge invariant evidence of the 
Abelian monopole condensation is still lacking. 

The aim of the present paper is to investigate the Abelian monopole condensation in 
pure lattice gauge SU(2) and SU(3) theories in a gauge-invariant way |12|]. To do this we 
introduce a disorder parameter defined in terms of a gauge-invariant thermal partition 
functional in presence of an external background field. 

The plan of the paper is as follows. In Section II we introduce the thermal partition 
functional, built up using the lattice Schrodinger functional []I3"|, |14T1 . In Section III we 
study the Abelian monopole condensation for finite temperature SU(2) lattice gauge the- 
ory. Section IV is devoted to the case of SU(3) gauge theory at finite temperature, where, 
according to the choice of the Abelian subgroup, different kinds of Abelian monopoles can 
be defined. Our conclusions are drawn in Section V. 



II. THE THERMAL PARTITION FUNCTIONAL 

To investigate the dynamics of the vacuum at zero temperature we introduced [T^,[T^| the 
gauge-invariant effective action for external static (i.e. time-independent) background 
field defined by means of the lattice Schrodinger functional: 

Z[U cxt ] = [vU e~ Sw , (2.1) 



where Sw is the standard Wilson action. The functional integration is extended over links 
on a lattice with the hypertorus geometry and satisfying the constraints 



U,(x)\ X4=0 = U; x \x). (2.2) 

In Equations (|2.1| ) and ( [2.2|) U^ xt (x) is the lattice version of the external continuum gauge 
field A ext (x) = A ext (x)X a /2: 



UJ? (x ) = P exp {iagj 1 dt A^x + tfi) ^ } 



(2.3) 




where P is the path-ordering operator and g the gauge coupling constant. 

The lattice effective action for the external static background field A ext (x) is given by 



T[A^] = --\n , (2.4) 



where L 4 is the extension in Euclidean time and Z(0) is the lattice Schrodinger functional, 
Eq. ( |2.1| ), without the external background field (U° xt = 1). It can be shown \T5\ that 

in the continuum limit r[A ext ] is the vacuum energy in presence of the background field 
A ext (x). 

We want now to extend our definition of lattice effective action to gauge systems at 
finite temperature. In this case the relevant quantity is the thermal partition function. 
In the continuum we have: 

Tr [e~P TH ] = J VA{A\e-^ H V\ A) , (2.5) 
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where (3? is the inverse of the physical temperature, H is the Hamiltonian, and V projects 
onto the physical states. As is well known, the thermal partition function can be written 

m 



as 



i J 4 M (/3 T ,2 f )=A M (0,S) 



On the lattice we have: 



Tr[e- fhH ]= f VU^x 4 ,x)e~ Sw . (2.7) 

Comparing Eq. Q2.7|) with Eqs. (|2.1|) and ([2.21) , we get: 

Tr [e~^ H ] = J VU,{x) Z[U^x)\ , (2.8) 

where Z[Un(x)] is the Schrodinger functional Eq. ( |2.1| ) defined on a lattice with L 4 = /3 T , 
with "external" links U^{x) at x 4 = 0. 

We are interested in the thermal partition function in presence of a given static back- 
ground field A ext (x). In the continuum this can be obtained by splitting the gauge field 
into the background field A ext (x) and the fluctuating fields r](x). So that we could write 
formally for the thermal partition function ^[A ext ]: 

Z T [A ext } = J Vff{A ext ,ff\e~ pTH V\A ext ,ff) . (2.9) 

The lattice implementation of Eq. fl2.9|) can be obtained from Eq. (|2.7p if we write 

U k (j3r, x) = £4(0, x) = U^(x)U k (x) , (2.10) 
where t/| xt (x) is given by Eq. ( |2.3|) and the f/fc(^)'s are the fluctuating links. Thus we get 

Z T [A cxt ] = J VU k (x)VU 4 (x)Z[U^(x),U k (x)], (2.11) 

where we integrate over the fluctuating links U k (x), while the U^ xt links are fixed. Note 
that in Eq. ( |2.11| ) only the spatial links belonging to the hyperplane x^ = are written 
as the product of the external link U% xt (x) and the fluctuating links U k (x). The temporal 
links Ui(x4 = 0,x) are left freely fluctuating. It follows that the temporal links U±{x) 
satisfy the usual periodic boundary conditions. We stress that the periodic boundary 
conditions in the temporal direction are crucial to retain the physical interpretation that 
the functional ^[y4 cxt ] is a thermal partition function. In the following the spatial links 
belonging to the time-slice 24 = will be called "frozen links" , while the remainder will 
be the "dynamical links" . 

From the physical point of view we are considering the gauge system at finite temperature 
in interaction with a fixed external background field. As a consequence, in the Wilson 
action Sw we keep only the plaquettes built up with the dynamical links or with dynamical 
and frozen links. With these limitations it is easy to see that in Eq. ( 2.11| ) we have 



Z 



Ul x \x),U k {x) =Z[Ul x \x)] . (2.12) 
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Indeed, let us consider an arbitrary frozen link U^. xt (x)Uk{x). This link enters in the 
modified Wilson action by means of the plaquette: 

P kA {x A = 0,x)= Tr tuF(g)U k fflUi(0,£+k)Ul{l,$Ut{0,$\ . (2.13) 

Now we observe that the link U^(0, x + k) in Eq. ( |2.13| ) is a dynamical one, i.e. we are 
integrating over it. So that, by using the invariance of the Haar measure we obtain 

P fc4 (x 4 = 0,f) =Tr jf4 cxt (f)£/4(0,f +k)Ui(l,x)Ul(0,x)} . (2.14) 



It is evident that Eq. (|2.14 ) in turns implies Eq. Q2.12 ). Then, we see that in Eq. ( [2.11 



the integration over the fluctuating links U(x) gives an irrelevant multiplicative constant. 
So that we have: 



A 



ext 



VUe~ Sw , (2.15) 

U k {p T ,x)=U k {0,x)=Ul^{x) 



where the integrations are over the dynamical links with periodic boundary conditions in 
the time direction. As concerns the boundary conditions at the spatial boundaries, we 
keep the fixed boundary conditions Uk(x, 24) = U^ dl {x) used in the Schrodinger functional 
Eq. (|2.1|) . Thus we see that, if we send the physical temperature to zero, then the thermal 
functional Eq. ( |2.15| ) reduces to the zero-temperature Schrodinger functional Eq. (|2.1| ) 



with the constraints Uk(x)\ X4=0 = L^ xt (x) instead of Eq. (|2.2|) . In our previous study jJOjl 



we checked that in the thermodynamic limit both conditions agree as concerns the zero- 
temperature effective action Eq. (|2.4j) . 



III. ABELIAN MONOPOLE CONDENSATION: SU(2) 

Let us consider the SU(2) pure gauge theory at finite temperature. We are interested in 
the thermal partition function Eq. ( p. 15 ) in presence of an Abelian monopole field. In the 



case of SU(2) gauge theory the maximal Abelian group is an Abelian U(l) group. Thus, 
in the continuum the Abelian monopole field turns out to be: 

gb a ( X ) = 5 a >* ^— . 3.1 

2 pKp| — x ■ n) 

where n is the direction of the Dirac string and, according to the Dirac quantization 
condition, integer. The lattice links corresponding to the Abelian monopole 

field Eq. ( |3.ip can be readily obtained as: 

UT\x) = Pexp S^ig J dt^b a k (x + tx k )^ , (3.2) 

where the cx a 's are the Pauli matrices. By choosing n = x 3 we get: 

U$(x) = cos[e lj2 (x)\ + ta 3 sin[0i, 2 (f)] , 

ut x Hx) = 1, 1 ' ' 
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with 



e 2 (x) = -+ 



"mon (x 2 ~ X 2 ) 1 

4 |*^mon| |-^mon| (^3 -^s) /q *\ 

n m0 n(xi-Xi) 1 



I ^mon | | ^mon | (^3 -^-3 ) 



In Equation fl3.4|) (X 1; X 2 , X 3 ) are the monopole coordinates and x mon — (x — X). In the 
numerical simulations we put the lattice Dirac monopole at the center of the time slice 
£4 = 0. To avoid the singularity due to the Dirac string we locate the monopole between 
two neighboring sites. We have checked that the numerical results are not too sensitive 
to the precise position of the magnetic monopole. 

According to the discussion in the previous Section we are interested in the thermal 
partition function 2r[A ext ] given by Eq. (|2.15| ). Note that we do not need to fix the gauge 



due to the gauge invariance of the thermal partition functional against gauge transforma- 
tions of the external background field. On the lattice the physical temperature T phys is 
given by 

Pr = L t , (3.5) 



T 

A phys 

where L t = L4 is the lattice linear extension in the time direction. In order to approximate 
the thermodynamic limit the spatial extension L s should satisfy 

L s ~>L t . (3.6) 

To this end we performed our numerical simulations on lattices such that 

^<4. (3.7) 

In the numerical simulations we impose periodic boundary conditions in the time direc- 
tion. As already discussed, at the spatial boundaries the links are fixed according to 
Eq. ( |3.3|) . This last condition corresponds to the requirement that the fluctuations over 
the background field vanish at infinity. 

Following the suggestion of Ref. |Tl]] we introduce the gauge-invariant disorder param- 
eter for confinement 

/x = e~ WW = , (3.8) 

where Zt[0] is the thermal partition function without monopole field (i.e. with n mon = 0). 

From Eq. ( |3.8| ) it is clear that F mon is the free energy to create an Abelian monopole. 
If there is monopole condensation, then F mon = and [1 = 1. To avoid the problem of 
measuring a partition function we focus on the derivative of the monopole free energy: 

d 

-^nion = ~QpFmon ■ (3-9) 

It is straightforward to see that -F^on ^ s gi ven by the difference between the average 
plaquette 



F' =V[<P >„ mon=0 - < P > nmo n ^o] , (3.10) 
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where V is the spatial volume. 

We use the over-relaxed heat-bath algorithm to update the gauge configurations. Sim- 
ulations have been performed by means of the APEIOO/Quadrics computer facility in Bari. 
Since we are measuring a local quantity such as the plaquette, a low statistics (from 2000 
up to 12000 configurations) is required in order to get a good estimation of F^ on . 

In Figure 1 we display the derivative of the monopole free energy versus j3 for n mon = 10 
on a lattice with L t = 4 and L s = 24. We see that -F mon vanishes at strong coupling and 
displays a rather sharp peak near (3 ±2 2.13. We expect that the peak corresponds to 
the finite temperature deconfinement transition. In Figure 1 we also display the absolute 
value of the Polyakov loop: 

x t=l 

and, indeed, we see that the peak corresponds to the rise of Polyakov loop. 
In the weak coupling region the plateau in F^ on indicates that the monopole free energy 
tends to the classical monopole action which behaves linearly in (3. To see this, we observe 
that deeply in the weak coupling region the lattice action should reduce to the classical 
action. In the naive continuum limit the classical action reads : 

1 rPr f 

Sciass = - J dx 4 J d 3 xB a (x) B a (x) (3.12) 

where B a (x) is the classical Abelian monopole magnetic field. Introducing an ultraviolet 
cutoff A = a/a, with a a constant and a the lattice spacing, and performing in Eq. (ft.12 ) 
the spatial integral over the volume V = L? s , we get: 

Sclass^ ^-^ + 0(1/^0). (3.13) 

" phys 

So that in the weak coupling region we have : 



From Figure 1 we see that Eq. ( ft.14 ) with a ±2 1.2 ( dashed line) describes quite well the 



numerical data in the relevant region. 

In order to determine the critical parameters and the order of the transition, we need 
to perform the finite size scaling analysis. We plan to do this in a future work. In this 
paper we restrict ourself to a preliminary qualitative analysis. In Figure 2 we compare 
the derivative of the monopole free energy on lattices with L t = 4 and L s = 24,48. We 
see that in the strong coupling -F mon agrees for the two lattices. On the other hand, in the 
weak coupling region the different values of the plateaus can be ascribed to finite volume 
effects. In the critical region we see that the peak increases. 

With the aim of obtaining the disorder parameter fi (Eq. (|3.8|) ) we perform the nu- 
merical integration of the monopole free energy derivative 

F mQ n(P)= I dp F' mon {(3') . (3.15) 
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In Figure 3 we show the disorder parameter fi versus (3 for lattices with L t = 4 and 
L s = 24,48. We see clearly that /i = 1 in the confined phase. In other words the 
monopoles condense in the vacuum. On the other hand, it seems that fi — > in the 
thermodynamic limit when j3 reaches the critical value . Indeed, by increasing the spatial 
volume of the lattice, the disorder parameter \x decreases faster toward zero. Moreover we 
see that the finite volume behavior of our disorder parameter is consistent with a second 
order deconfinement phase transition. 

It worthwhile to comment on the finite volume effects. As a matter of fact, it appears 
that, even though the spatial volume of our larger lattice looks enormous, we gain a 
rather small increase in the peak value of the monopole free energy derivative. This can 
be understood by observing that, due to our peculiar conditions at the spatial boundaries, 
the dynamical volume is smaller than the geometrical one. Moreover, it is well known 
that the fixed boundary conditions for the gauge fields lead to more severe finite volume 
effects with respect to the usual periodic boundary conditions. So that, to reach the 
thermodynamic limit we must simulate our gauge system on lattices with very large 
spatial volumes. We stress again that the precise determination of the critical parameters 
requires a finite size scaling which will be presented elsewhere. 

IV. ABELIAN MONOPOLE CONDENSATION: SU(3) 

In the case of SU(3) gauge theory, the maximal Abelian group is U(l)xU(l). Therefore 
we have two different types of Abelian monopole. Let us consider, firstly, the Abelian 
monopole field given by Eq. ( |3.1| ), which we call the T 3 Abelian monopole. The lattice 
links are given by 



U$(x) 



3»01,2(e) 








e -i0i j2 (x) Q 

1 



(4.1) 



with 9i t 2(x) defined in Eq. (|3.4| ). The second type of independent Abelian monopole can 
be obtained by considering the diagonal generator T§. In this case we have the T§ Abelian 
monopole: 



U$(x) 
















-2ie 1<2 (x) 



(4.2) 



with 



1 


^■mon 


(x 2 - X 2 ) 




1 




4 


I ^mon | 


| ^mon | 


- (x 3 - X 3 ) 


i 




(xi-Xx) 




1 


71 


4 


1 | 


I ^mon | 


- (x 3 - X 3 ) 



(4.3) 



Obviously, the lattice links Eq. ([4.2|) corresponds now to the continuum gauge field 



gb a (x) = 6' 



a a n mon 



x x n 



2 151(1x1 — x • n) 



(4.4) 
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The other Abelian monopoles can be generated by considering the linear combination of 
the T 3 and T 8 generators. In particular we have considered the T 3a Abelian monopole 
corresponding to the following linear combination []TT| of A 3 /2 and Ag/2: 



T - — + — — 

3a ~~ 2 2 2 2 




\ 

-I. 



(4.5) 



In Figure 4 we compare the free energy monopole derivative for the T 3 , T 3a and T 8 Abelian 
monopoles for the lattice with L t = 4 and L s = 32. We see that the T 3 and T 3a Abelian 
monopoles agree within statistical errors in the whole range of (3. On the other hand 
the Tg Abelian monopole displays a signal about a factor two higher in the peak region. 
This is at variance of previous studies 0] which find out that the disorder parameters 
for the three Abelian monopoles defined by means of the Polyakov projection coincide 
within statistical errors. This result is quite interesting, for it suggests that in the pattern 
of dynamical symmetry breaking due to the Abelian monopole condensation the color 
direction 8 is slightly preferred. 

Let us consider, now, in detail the T 8 Abelian monopole. In Figure 5 we report the 
derivative of the monopole free energy versus /3 for the lattice with L t = 4 and L s = 32. 
We also display the absolute value of the Polyakov loop: 



t=i 



We see that F^ on behaves like in the SU(2) case. Indeed, the free energy monopole 
derivative is zero within errors in the strong coupling region, while it display a sharp peak 
in correspondence of the rise of the Polyakov loop. In the weak coupling region -F^on * s 
almost constant. The value of the plateau correspond to the classical action Eq. ( |3.12| ) 
which in the present case gives: 

(X7T 3 

Sciass ^ < on + 0(1/L s a) , (4.7) 

J-^ ■*■ phys 



so that 



The dashed line in Figure 5 in the weak coupling region corresponds to Eq. Q4.8|) with 
a ~ 2.0. 

As in the SU(2) theory we find that by increasing the spatial volume the peak increases 
(see Figure 6). Our data do not show a measurable shift of the peak. We feel that this 
is a manifestation of the first order nature of the 577(3) deconfinement transition. This 
is confirmed if we look at the disorder parameter fi. In Figure 7 we show the disorder 
parameter fi versus Q for the L t = 4 and L s = 32,48 lattices. Again we see that the 
disorder parameter fi is different from zero in the confined phase and decreases towards 
zero in the thermodynamic limit when we approach the critical coupling. Moreover our 
numerical results suggest that by increasing the spatial volume the two curves cross. This 
is precisely the finite volume behavior expected for the order parameter in the case of a 



first order phase transition [18 



S 



V. CONCLUSIONS 



In this paper we have investigated the Abelian monopole condensation in the finite tem- 
perature SU(2) and SU(3) lattice gauge theories. By means of the lattice thermal partition 
functional we introduce a disorder parameter which signals the Abelian monopole conden- 
sation in the confined phase. By construction our definition of the disorder parameter is 
gauge invariant, so that we do not need to perform the Abelian projection. Our numerical 
results suggest that the disorder parameter \x is different from zero in the confined phase 
and tends to zero when approaching the critical coupling in the thermodynamic limit. We 
point out that in our approach the precise determination of the critical parameters could 
be obtained by means of a finite size scaling analysis. However, our results are consistent 
with a second order deconfining phase transition in the case of the SU (2) gauge theory. 
On the other hand, in the case of SU (3) the disorder parameter \x displays the finite-size 
behavior expected for a first order transition. It is clear that the finite size analysis in 
the critical region requires a separate study with both better statistic and larger lattice 
volumes. Remarkably, in the case of SU(3) gauge theory, where there are two independent 
Abelian monopole fields related to the two diagonal generators of the gauge algebra, we 
find that the non perturbative vacuum reacts moderately strongly in the case of the Tg 
Abelian monopole. We feel that this last result should be useful in the theoretical efforts 
to understand the pattern of symmetry breaking in the deconfined phase of QCD. 
In conclusion we stress that our approach, while keeping the gauge invariance, can be 
readily extended to incorporate the dynamical fermions. We hope to present results in 
this direction in a future study. 
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Figure 1: The derivative of the SU(2) monopole free energy versus j3 (Eq. ( p. 10 )) for 
n mon = 10 on a lattice with L t = 4 and L s = 24 (open circles), with the absolute value of 
the Polyakov loop (full circles). The dashed line is Eq. (|3.13|) . 
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Figure 2: The derivative of the SU(2) monopole free energy versus (3, (Eq. (|3.10|) ), on 
lattices with L t = 4 and L s = 24 (open circles), and L s = 48 (full circles). 
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Figure 3: The disorder parameter /x (Eq. ( |3.8| )) for the SU(2) monopoles versus (3 for 
lattices with L t — 4 and L s = 24 (open circles), or L s = 48 (full circles). 
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Figure 4: The derivative of the SU(3) monopole free energy in the case of T3 (circles), T^ a 
(squares), and T 8 (diamonds) Abelian monopoles. 
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Figure 5: The derivative of the SU(3) monopole free energy for the T 8 Abelian monopole 
(open circles) versus (3 with the absolute value of the Polyakov loop (full circles). The 
dashed line is Eq. Q4.8Q . 
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Figure 6: The derivative of the T 8 Abelian monopole free energy versus (3, on lattices with 
L t — 4 and L s = 24 (open circles), and L s = 48 (full circles). 
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Figure 7: The disorder parameter \i (Eq. ( |3.8|) ) for the SU(3) T 8 Abelian monopoles versus 
f3 for lattices with L t = 4 and L s = 24 (open circles), and L s = 48 (full circles). 
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